Open Access
How to translate text using browser tools
30 January 2001 A Molecular Perspective on the Phylogeny of the Girdled Lizards (Cordylidae, Squamata)
DARREL FROST, DANIEL JANIES, P. LE FRAS N. MOUTON, TOM TITUS
Author Affiliations +
Abstract

Mitochondrial DNA sequences were obtained for 16 species representing all nominal genera of Cordylidae (Platysaurus, Chamaesaura, Cordylus, and Pseudocordylus). Gerrhosauridae and Teiidae were used as first and second outgroups. Results indicate that the oviparous Platysaurus is the sister taxon of the remaining cordylids (all of which are ovoviviparous). Within the ovoviviparous group Cordylus is paraphyletic with respect to Chamaesaura and Pseudocordylus. No evidence of Pseudocordylus monophyly was discovered. The three species of Chamaesaura and the seven species of Pseudocordylus are transferred to Cordylus to render a monophyletic taxonomy.

INTRODUCTION

The Cordylidae, girdled lizards, form a monophyletic group composed of more than 50 species arranged in four nominal genera (Loveridge, 1944; Branch, 1998): Platysaurus, a group of 15 species most notable for their extremely flattened shape; Pseudocordylus, a group of seven species most notable for their apparent intermediacy in flatness between Platysaurus and deeper-bodied Cordylus; Cordylus, a group of 31 species of diverse morphology from extremely spiny to relatively smooth scaled; and finally, Chamaesaura, a group of three serpentiform species with reduced limbs. The Cordylidae are restricted to the southern subcontinent of Africa, except for a few species that extend into central and East Africa.

Knowledge of cordylid phylogeny is in its infancy, even though its near outgroups and its monophyly have been known for a considerable length of time (McDowell and Bogert, 1954; Estes et al., 1988). Most of the systematic work on cordylid lizards has been focused primarily on species limits and geographic variation (e.g., FitzSimons, 1943; Loveridge, 1944), although recent work has also attempted to place species within a phylogenetic context (e.g., Mouton, 1990; Mouton and Van Wyk, 1989, 1994, 1995; Harvey and Gutberlet, 1995). The most comprehensive attempt at phylogeny reconstruction within the Cordylidae is that of Lang (1991), and although that paper provides considerable information on morphological variation within the Cordylidae and its first taxonomic outgroup, Gerrhosauridae, it only provides a beginning to our understanding of phylogenetics within the group. Genera of unsubstantiated monophyly were used as terminal taxa in Lang's revision, even though evidence to the contrary existed then in the form of considerable interspecific variation in many characters showing coherent distributions across generic boundaries. Unfortunately, a data matrix for species was not provided even though the apomorphy and change lists provided in appendices suggest that such a matrix exists. Many of the characters discussed by Lang provide tantalizing evidence for the paraphyly of such nominal genera as Cordylus and Pseudocordylus, even though his evidence, reified through the assumption of generic monophyly, necessarily appears to support monophyly of those taxa.

Our initial intentions in this project were (1) to use the morphological data provided by Lang (1991) and Harvey and Gutberlet (1995) as external criteria to evaluate differences in phylogenetic results due to sensitivity to alignment parameters (i.e., to search for topologies that minimize incongruence among data sets due to alignment parameters—sensu Wheeler, 1995) and (2) to perform a character congruence study to enhance our understanding of cordylid relationships. With respect to Lang's (1991) study, this has proven impossible for reasons of incompatible taxon sampling as well as our inability to extract species-level information from that revision with confidence. Coupled with our doubt regarding the monophyly of at least two of the four genera (i.e., Pseudocordylus and Cordylus) assumed to be monophyletic by Lang, these difficulties have proven insurmountable. With respect to the scale surface ultrastructure study of Harvey and Gutberlet (1995) there were very few characters addressed (and, to foreshadow our results, only three characters in conflict with the molecular evidence). The third author of this paper, P. Mouton, has for some time been engaged in a survey of the taxonomy, morphology, and anatomy of the group (e.g., Mouton, 1985, Mouton et al., 1992; Brody et al., 1993; Mouton and Van Wyk, 1994, 1997). A more inclusive publication of the morphological data will be provided by Mouton at a later date. Therefore, the purpose of this contribution is not to report on the morphology of the cordylids, but to report on an analysis of DNA sequences with comments on what these data suggest with respect to reported morphological diversification within the Cordylidae. We will mention the morphological evidence of Lang (1991) and Harvey and Gutberlet (1995), but only within a narrative context.

METHODS

Sampling

Mitochondrial DNA sequences were obtained for 16 terminals representing all nominal genera as well as the putative first and second functional outgroup. (The first taxonomic outgroup, the Gerrhosauridae, is represented by Gerrhosaurus typicus and Tetradactylus seps; the second taxonomic outgroup is represented by a teiid, Cnemidophorus sexlineatus.) See appendix 1 for voucher numbers and associated GenBank accession numbers.

DNA Amplification and Sequencing

Mitochondrial genes encoding the 12S rDNA, valine tDNA, and the 5′ end of the 16S rDNA were amplified using the polymerase chain reaction. Genomic DNA extraction, primers, and amplification protocols were identical to those in Titus and Frost (1996). Amplified DNA was electrophoresed on 1% agarose gels and purified for sequencing using the Geneclean II kit (Bio 101, Inc.). Thermal cycle sequencing was done following Titus and Frost (1996). Automated sequencing was performed in the University of Oregon Molecular Biology Sequencing Facility utilizing the Big Dye Terminator Cycle Sequencing Kit with AmpliTaq FS (Perkin-Elmer) and an ABI PRISM 377 DNA Sequencer (Perkin-Elmer) following manufacturer's specifications.

Alignment and Analysis

Sequences were not aligned prior to a parsimony analysis, which is the usual procedure for most sequence analyses. Instead, in this case we employed direct optimization (Wheeler, 1996) as implemented by the computer program POY (Gladstein and Wheeler, 1996–1999). (For summaries of this method and access to the program and command scripts see ftp.amnh.org.) Although there is obviously considerable structure in the data as suggested in Results and Discussion by the Bremer values, the g1 statistic as a significance test was not employed for reasons detailed by Kallersjö et al. (1992). A maximum likelihood approach has also not been employed because these have been demonstrated to converge on parsimony estimates when there is no assumption of a common model of evolution across sites (Tuffley and Steel, 1997).

Gap, transversion, and transition costs to the analysis procedure were applied in ratios of 1:2:1 and 1:1:1, for reasons detailed by Titus and Frost (1996) and to provide a conservative estimate of cladistic structure. Although we reject the notion that support statistics are measures of accuracy, there is no doubt that they are quite useful as practical indicators of whether we are making progress toward efficient and consistent description. For this reason, Bremer indices (Bremer, 1994) are provided.

RESULTS AND DISCUSSION

Topologies for gap, transversion, and transition costs of 1:1:1 and 1:2:1 are provided in figures 1 and 2, respectively. Length for the two trees was, respectively, 1501 and 2257, with consistency and retention indices of, respectively, 0.22 and 0.80, and 0.25 and 0.80. Bremer (= decay) indices (Bremer, 1994) are provided with the relevant stems on figures 1 and 2. A strict consensus of these two topologies is provided in figure 3. The genus-terminal tree of cordylids provided by Lang (1991, his fig. 17) and a summary of our results are provided in figure 4 for comparison.

Cost parameters showed some effect with a difference of 2× in the weighting of transversions over transitions and gaps. Only one stem with a Bremer index of 10 was affected by this change (affecting the relative position of Cordylus coeruleopunctatus), and monophyly of the putative sister-taxon, Gerrhosauridae, is not supported by the molecular data. This is a surprising result and may have been caused by our limited sampling. However, the Bremer indices are high enough that reevaluation of the previously published evidence for gerrhosaurid monophyly would be prudent.

Monophyly of the Cordylidae sensu Lang (1991), always noncontroversial, is affirmed. Lang (1991) summarized a number of morphological synapomorphies for the group and it must be one of the most highly corroborated taxa in lizard systematics.

Only one species of Chamaesaura is included in the analysis, but the extraordinary morphology of this taxon strongly suggests monophyly of the three included species. The suggestion of the molecular data—that Chamaesaura rests within nominal Cordylus—is a novel hypothesis and at variance with the previous suggestion by Lang (1991) that Chamaesaura is the sister taxon of other cordylids. Beyond a large number of generic apomorphies, Lang (1991) suggested five characters for which Chamaesaura was plesiomorphic with respect to the remaining cordylids, but he also mentioned several other features of Chamaesaura that were found in some Cordylus (e.g., his S1-position of nasal scales, S3-condition of postnasal scale). Further, Chamaesaura is so highly apomorphic, correlated with its serpentiform lifestyle, that our molecular data suggest that careful evaluation of the morphological data is required. Note, however, that the molecular data placing Chamaesaura within “Cordylus” are corroborated by at least one anatomical-ecophysiological characteristic. Mouton and Van Wyk (1997) noted that outgroup comparison would judge ovoviviparity within Chamaesaura, Cordylus, and Pseudocordylus as a synapomorphy and that oviparity in Platysaurus (the only ovoviparous cordylid) would be counted as a plesiomorphy, a hypothesis congruent with our molecular results.

Although we have only two representatives of Platysaurus (P. monotropis and P. intermedius rhodesianus) we have no reason to doubt its monophyly. In Lang's (1991) cladogram, Platysaurus is diagnosed by seven apomorphies. The molecular data place Platysaurus as the sister-taxon of the remaining cordylids, including Chamaesaura. This, like the placement of Chamaesaura within Cordylus based on molecular evidence, is surprising inasmuch as what appears to be a morphocline—from the smooth-scaled, flattened crevice-dwellers (Platysaurus) through the somewhat less flattened Pseudocordylus to variably spiny Cordylus—had initially appeared to us as a phylogenetic transition (with the Chamaesaura morphotype independently derived from the Cordylus morphotype). The molecular data, like the ecophysiological data of Van Wyk (1989) and Mouton and Van Wyk (1997), suggest that this is not the case, although this placement within Cordylus resolves the problem of reversal from ovoviviparity (Chamaesaura, Cordylus, and Pseudocordylus) to oviparity (Platysaurus and cordylid outgroups) required by Lang's placement.

No evidence of monophyly of Pseudocordylus is provided by the molecular data. Two clades of Pseudocordylus are recognized by our analysis, the first composed of Ps. microlepidotus and Ps. melanotus, and the second composed of Ps. capensis and Ps. nebulosus. Regardless of their similarity, there is no molecular support for a special relationship of Pseudocordylus with Platysaurus and considerable evidence against it. Lang (1991) supported the monophyly of Pseudocordylus with three characters: (a) nasal scales in contact (his S1.1; even though he suggested that Platysaurus, the putative sister taxon, was highly variable in this characteristic), (b) anterior parietal scales in contact (his S11.1; although he suggested that Platysaurus, the putative sister taxon in his analysis, showed contact in most species and variability in Cordylus was high), and (c) lingual exposure of the angular (his C24.0, although he noted on p. 185 that it is variable among unspecified species of Pseudocordylus). In summary, we view the evidence for monophyly of Pseudocordylus to be equivocal and consider Pseudocordylus a candidate for being polyphyletic. Further, at least in this study, evidence tying either group of Pseudocordylus to Platysaurus appears to be lacking.

No evidence of monophyly of Cordylus is suggested by the molecular evidence. Indeed, there is very strong evidence of paraphyly with respect to Pseudocordylus and Chamaesaura. That Cordylus is paraphyletic with respect to Pseudocordylus is not surprising. Lang (1991) suggested five synapomorphies of Cordylus: S21.1 (keeled scales on bottom of foot; not scorable in Chamaesaura and variable in Platysaurus and Pseudocordylus), S24.1 (keeled subdigital lamella; not scorable in Chamaesaura and variable in Pseudocordylus), C28.0 (posteriorly tapered retroarticular process [a character influenced by degree of skull flattening]), P1.2 (intercentra fused to posterior vertebrae [sutured in Chamaesaura]; and M3.0 (osteoderm distribution, osteoderms apparently always reduced in flattened crevice dwellers). Nevertheless, even though Cordylus was assumed to be monophyletic on the basis of these characters (and for purposes of his intergeneric analysis), Lang (1991) noted at least nine features that are variable among species of Cordylus and that extend to other cordylid taxa.

Differences due to transversion weighting are evident, particularly in the placement of the Cordylus vittifer group, C. tropidosternum, C. warreni, C. coeruleopunctatus, and the relative placement of Chamaesaura within “Cordylus” clade.

Clades that remain the same among our analyses within nominal Cordylus are: (a) a clade of Cordylus oelofseni, C. cordylus, C. niger, C. macropholis, and C. minor; (b) a Cordylus vittifer, C. peersi, C. lawrenci clade; (c) a clade composed of C. jordani and C. polyzonus. These associations are hardly surprising. Mouton (1990) first noted the association of C. oelofseni, C. cordylus, and C. niger. Loveridge (1944) considered these three species to be parts of a single larger species, C. cordylus, but he had also included a number of other populations as well, including C. lawrenci, judged to be distant from this group on the basis of our data. The close association of C. jordani and C. polyzonus has also been long recognized (Loveridge, 1944). Nevertheless, the scale ultrastructure data presented by Harvey and Gutberlet (1995), and discussed below, suggest that continued careful evaluation of relationships is warranted.

The scale surface ultrastructure data of cordylids were cast by Harvey and Gutberlet (1995) into 10 characters of transformation and mapped upon the cladogram of Lang (1991), although Harvey and Gutberlet discussed some additional refinement within Cordylus. When the accuracy of Lang's (1991) generic arrangement is not assumed, very little discordance with Harvey and Gutberlet's data and our results is discovered and what is found has no impact on our hypothesis of Cordylus paraphyly with respect to Pseudocordylus and Chamaesaura. Indeed, when their data are analyzed independently they show considerable concordance with our molecular results (Stems A and B, fig. 5; strict consensus of a sample of 16383 trees; length = 11, CI = 0.9091; without uninformative characters 0.8750). What discordance exists is solely in regard to the relative placement of various species of Cordylus (fig. 5, stems C and D). The presence on dorsal scales of the macrohoneycomb oberhautchen (character #1 of Harvey and Gutberlet), thought to be lost in Platysaurus when mapped on Lang's arrangement is plesiomorphically absent in our cladogram in Platysaurus and its presence is a synapomorphy of the remaining cordylids. (In other words, this character corresponds in placement with the transition from oviparity to ovoviviparity.) Where their evidence does show discordance with the molecular evidence is in their character #6 (flaplike free margins associated with the cell-ridge system; fig. 5, Stem C), and their character #7 (ridges in the center of oberhautchen cells in C. polyzonus, C. niger, and C. minor; fig. 5, stem D). Although discordance between ultrastructural character evidence and the molecular evidence within Cordylus bears further scrutiny, in this case the molecular evidence is rather strong, suggesting the existence of more homoplasy in scale surface ultrastructure than previously suspected.

Clearly a study of character congruence of both morphology and molecules of species (not of higher taxa as was done by Lang, 1991) would provide a stronger basis on which to optimize transversion, transition, and gap cost parameters and would further elucidate the relative relationships within cordylids. Nevertheless, on the basis of the molecular evidence alone, it is evident that “Cordylus” is paraphyletic with respect to Chamaesaura and Pseudocordylus, that Pseudocordylus is dubiously monophyletic, and that the oviparous Platysaurus is the sister taxon of the other, ovoviviparous, cordylids. On this basis we regard Pseudocordylus and Chamaesaura to be junior synonyms of Cordylus. This affects only 10 species names and formulates a conservative taxonomy consistent with the state of knowledge of cordylid relationships.

Acknowledgments

This contribution was supported by NSF Grant DEB-9220870 to D. Frost and T. Titus, by NASA grant NAG5–8443 to D. Frost and W. Wheeler, and by an AMNH postdoctoral fellowship to D. Janies. We especially appreciate thoughtful criticism from A. Bauer, J. A. Campbell, J. Faivovich, T. Grant, M. Harvey, J. Losos, and W. Wheeler.

REFERENCES

1.

B. Branch 1998. Field guide to snakes and other reptiles of southern Africa. Sanibel: Ralph Curtis Books. Google Scholar

2.

K. Bremer 1994. Branch support and tree stability. Cladistics 10:295–304. Google Scholar

3.

L. S. Brody, Ple F N. Mouton, and W. S. Grant . 1993. Electrophoretic analysis of the Cordylus cordylus species group in the southwestern Cape, South Africa. Amphibia-Reptilia 14:19–33. Google Scholar

4.

R. Estes, K. de Queiroz, and J. Gauthier . 1988. Phylogenetic relationships of the Squamata. In R. Estes and G. Pregill (eds.), Phylogenetic relationships of the lizard families: 119–291. Stanford: Stanford Univ. Press. Google Scholar

5.

V. FitzSimons 1943. The lizards of South Africa. Pretoria. Google Scholar

6.

D. Gladstein and W. Wheeler . 1996–1999. POY. Software available at ftp.amnh.org in directory pub/people/wheeler/poy.  Google Scholar

7.

M. B. Harvey and R. L. Gutberlet, Jr . 1995. Microstructure, evolution, and ontogeny of scale surfaces in cordylid and gerrhosaurid lizards. J. Morphol 226:121–139. Google Scholar

8.

M. Kallersjö, J. S. Farris, A. G. Kluge, and C. Bult . 1992. Skewness and permutation. Cladistics 8:275–287. Google Scholar

9.

M. Lang 1991. Generic relationships within Cordyliformes (Reptilia: Squamata). Bull. Inst. R. Sci. Nat. Belg 61:121–188. Google Scholar

10.

A. Loveridge 1944. Revision of the African lizards of the family Cordylidae. Bull. Mus. Comp. Zool 95:1–118. 12 pl. Google Scholar

11.

S. B. McDowell and C. M. Bogert . 1954. The systematic position of Lanthonotus and the affinities of the anguinomorphan lizards. Bull. Am. Mus. Nat. Hist 105:1–142. Google Scholar

12.

Ple F N. Mouton 1985. Comments on the Cordylus cordylus complex in the southwestern Cape. J. Herpetol. Assoc. Africa 31:16–18. Google Scholar

13.

Ple F N. Mouton 1990. Taxonomic status of the melanistic forms of the Cordylus cordylus complex (Reptilia: Cordylidae) in the south-western Cape, South Africa. S. Afr. J. Zool 25:31–38. Google Scholar

14.

Ple F N. Mouton and J. H. Van Wyk . 1989. Cordylus minor: A valid species of South African lizard (Reptilia: Cordylidae). S. Afr. J. Zool 24:322–328. Google Scholar

15.

Ple F N. Mouton 1994. Taxonomic status of geographical isolates in the Cordylus minor complex (Reptilia: Cordylidae): A description of three new species. J. Herpetol. Assoc. Africa 43:6–18. Google Scholar

16.

Ple F N. Mouton 1995. A new crag lizard from the Cape Folded Mountains in South Africa. Amphibia-Reptilia 16:389–399. Google Scholar

17.

Ple F N. Mouton 1997. Adaptive radiation in cordyliforms lizards: an overview. Afr. J. Herpetol 46::76–88. Google Scholar

18.

Ple F N. Mouton, J. D. Visser, and J. H. Van Wyk . 1992. The status of the races of the graceful crag lizard, Pseudocordylus capensis, from South Africa. Amphibia-Reptilia 13:109–119. Google Scholar

19.

T. Titus and D. R. Frost . 1996. Molecular homology assessment and phylogeny in the lizard family Opluridae (Squamata: Iguania). Mol. Phylogenet. Evol 6:49–62. Google Scholar

20.

C. Tuffley and M. Steel . 1997. Links between maximum likelihood and maximum parsimony under a simple model of site substitution. Bull. Math. Biol 59:581–607. Google Scholar

21.

J. H. Van Wyk 1989. The female reproductive cycle of the lizard, Cordylus polyzonus polyzonus (Sauria: Cordylidae) in the Orange Free State. S. Afr. J. Zool 24:263–269. Google Scholar

22.

W. C. Wheeler 1995. Sequence alignment, parameter sensitivity, and the phylogenetic analysis of molecular data. Syst. Biol 44:321–332. Google Scholar

23.

W. C. Wheeler 1996. Optimization alignment: The end of multiple sequence alignment in phylogenetics? Cladistics 12:1–10. Google Scholar

Appendices

APPENDIX 1: SPECIMEN VOUCHERS

All specimens deposited in the Ellermann Museum, University of Stellenbosch, Stellenbosch, Rep. South Africa, except for the single specimen of Cnemidophorus sexlineatus that is unnumbered in the Department of Herpetology, University of Kansas. Cordylidae: Cordylus anguinus (formerlyChamaesaura anguina) (USEC-3503; GenBank BankIt321210-AF236042); C. capensis (formerly Pseudocordylus capensis) (USEC-H2486; GenBank BankIt319965-AF236022); C. cataphractus (USEC-H3490; GenBank BankIt320915-AF236034); C. coeruleopunctatus (USEC-H3497; GenBank BankIt321214-AF236043); C. cordylus (USEC-H3488; GenBank BankIt320904-AF236027); C. giganteus (USEC-H3494; GenBank BankIt321163-AF236036); C. jordani (USEC-H3492; GenBank BankIt320897-AF236024); C. lawrenci (USEC-H3506; GenBank BankIt321218-AF236046); C. macropholis (USEC-H3500; GenBank BankIt321217-AF236045); C. melanotus (formerly Pseudocordylus melanotus) (USEC-H3483; GenBank BankIt320905-AF236028); C. m. microlepidotus (formerly Pseudocordylus m. microlepidotus) (USEC-H3484; GenBank BankIt320910-AF236029); C. m. namaquensis (formerly Pseudocordylus m. namaquensis) (USEC-H3485; GenBank BankIt321208-AF236041); C. minor (USEC-H3498; GenBank BankIt320912-AF236032); C. niger (USEC-3489; GenBank BankIt320899-AF236025); C. nebulosus (formerly Pseudocordylus nebulosus) (USEC-H3487; GenBank BankIt320913-AF236033); C. oelofseni (USEC-H3499; GenBank BankIt321215-AF236044); C. peersi (USEC-H3495; GenBank BankIt320911-AF236031); C. polyzonus (USEC-H3491; GenBank BankIt320908-AF236030); C. tropidosternum (USEC-H3496; GenBank BankIt320859 AF236023); C. vittifer (USEC-H3507; GenBank BankIt320902-AF236026); C. warreni (USEC-H3493; GenBank BankIt320914-AF236035); Platysaurus intermedius rhodesianus (USEC-H3501; GenBank BankIt321165-AF236037); P. monotropis (USEC-H3502; GenBank BankIt321169-AF236038); Gerrhosauridae: Gerrhosaurs typicus (USEC-H3504; GenBank BankIt321222-AF236039); Tetradactylus seps (USEC-H3505; GenBank BankIt321207-AF236040). Teiidae: Cnemidophorus sexlineatus (BankIt321219-AF236047).

 Fig. 1. 

Tree based on relative cost parameters (transitions: transversions: gaps) at 1: 1: 1. Numbers on internal stems are Bremer values

i0003-0082-3310-1-1-f01.gif

 Fig. 2. 

Tree based on relative cost parameters (transitions: transversions: gaps) at 1: 2: 1. Numbers on internal stems are Bremer values

i0003-0082-3310-1-1-f02.gif

 Fig. 3. 

Strict consensus of trees shown in figures 1 and 2.

i0003-0082-3310-1-1-f03.gif

 Fig. 4. 

Summary trees of relationships suggested by Lang (1991) and here

i0003-0082-3310-1-1-f04.gif

 Fig. 5. 

Tree based on scale structure evidence of Harvey and Gutberlet (1995). A: macrohoneycomb present on venter; B: macrohoneycomb present on dorsal scales; C: flaplike free margins associated with the cell-ridge system; D: short ridgelike projections in the center of the oberhautchen cells

i0003-0082-3310-1-1-f05.gif
DARREL FROST, DANIEL JANIES, P. LE FRAS N. MOUTON, and TOM TITUS "A Molecular Perspective on the Phylogeny of the Girdled Lizards (Cordylidae, Squamata)," American Museum Novitates 2001(3310), 1-12, (30 January 2001). https://doi.org/10.1206/0003-0082(2001)310<0001:AMPOTP>2.0.CO;2
Published: 30 January 2001
Back to Top